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ABSTRACT 

Experimental data and numerical simulations of 
low-pressure turbines have shown that unsteady 
blade row interactions and separation can have a sig- 
nificant impact on the turbine efficiency. Measured 
turbine efficiencies at takeoff can be as much as two 
points higher than those at cruise conditions. Several 
recent studies have revealed that the performance of 
low-pressure turbine blades is a strong function of 
the Reynolds number. In the current investigation, 
experiments and simulations have been performed to 
study the behavior of a low-pressure turbine blade at 
several Reynolds numbers. Both the predicted and 
experimental results indicate increased cascade losses 
as the Reynolds number is reduced to the values as- 
sociated with aircraft cruise conditions. In addition, 
both sets of data show that tripping the boundary 
layer helps reduce the losses at lower Reynolds num- 
bers. Overall, the predicted aerodynamic and perfor- 
mance results exhibit fair agreement with experimen- 
tal data. 

NOMENCLATURE 

b x Axial chord 
Cj Skin friction coefficient 
Cp (P n -P)/(l/2 P U\) 

M Mach number 
P Pressure 
P.S. Pressure surface 
P t Total Pressure 


‘Associate Professor, Senior Member AIAA. 
a Research Engineer, Member ATAA. 
b Associate Professor, Senior Member ATAA. 
c Aerospace Engineer, Senior Member AIAA. 
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Re 

Reynolds number (axial chord, inlet vel 

S 

Entropy, arc-length distance 

s.s. 

Suction surface 

u 

Free stream velocity 

X 

Axial-direction 

Q 

Absolute reference frame flow angle 

8* 

Displacement thickness 

Tu 

Turbulence intensity 

cj 

Loss, (Pti - P t2 )/{l/2pUi 2 ) 


SUBSCRIPTS 

expt 

Experiment 

pred 

Predicted 

1 

Blade inlet 

2 

Blade exit 


INTRODUCTION 


Experimental data from jet-engine tests have in- 
dicated that unsteady blade row (wake) interactions 
and separation can have a significant impact on the 
efficiency of turbine stages. The effects of these in- 
teractions can be intensified in low-pressure turbine 
stages because of the low Reynolds number operat ing 
environment. Measured turbine efficiencies at takeofT 
can be as much as two points higher than those at 
cruise conditions [1]. Thus, during the last decade a 
significant amount of effort has been put into deter- 
mining the effects of transition and turbulence on the 
performance of low pressure turbine stages. Experi- 
mental investigations have been performed, for exam- 
ple, by Hodson et al. [2, 3], Halstead et al. [4, 5], Qiu 
et al. [6], Sohn et al. [7] and Boyle et al. [8]. These 
investigations have helped identify/clarify the roles 
that factors such as the Reynolds number, free stream 


turbulence intensity, pressure gradient and curvature 
have in the generation of losses. In particular, it has 
been determined that [1]: 

• At low to moderate Reynolds numbers there is 
a laminar region extending some distance from 
the leading edge. 

• The boundary layer may separate, particularly 
on the suction surface of the blade. Separation 
may occur in the form of a closed bubble, or 
as massive separation with no reattachment (re- 
sulting in large losses). The pressure surface may 
have cove separation, and small separation bub- 
bles may exist near t he leading edge. The separa- 
tion bubbles often originate in transitional flow, 
while reattachment is usually in turbulent flow. 

• The interaction of incoming wakes with the 
boundary layer often creates a convected tran- 
sitional or turbulent patch, which is trailed by a 
“calmed” region. The calmed region is a relax- 
ation region between the patch and the laminar 
boundary layer. 

More recently, Lake et ai [9, 10] performed experi- 
ments for the PAK-B blade, which was developed by 
Pratt & Whitney to operate at low speeds while at- 
taining loadings similar to those achieved at engine 
conditions. The study of Lake presented boundary 
profiles and loss data for three Reynolds numbers. In 
addition, the use of dimples and grooves for reducing 
losses at low Reynolds numbers was investigated. 

In parallel to the experimental investigations, there 
have been significant analytical efforts to improve the 
modeling of transition. Examples of such efforts in- 
clude the works of Mayle [11], Reshotko [12] and 
Gostelow et al. [13, 14]. These newer models show 
promise of providing accurate transition predictions 
over a wide rangle of flow conditions [15], although 
they have yet to be implemented into the numeri- 
cal flow analyses used by the turbine design commu- 
nity. Some recent computational investigations of in- 
terest include the works of Chernobrovkin and Lak- 
shminarayana [16], Kang and Lakshminarayana [17], 
Huang and Xiong [18] and Dorney and Ashpis [19]. 
The work of Dorney and Ashpis, similar to the ex- 
perimental efforts of Lake et ai , involved the PAK-B 
geometry (although at higher Mach numbers). 

The focus of the current effort has been to study 
the performance of low-pressure turbine airfoils by co- 
ordinating experiments and simulations for the PAK- 
B geometry. The PAK-B airfoil design was chosen 
because it generates a loading profile similar to the 
rotors found in modern aircraft engines. The com- 
parison of the experimental and numerical results in- 


cludes boundary layer profiles, airfoil loadings and 
losses. 

EXPERIMENTAL CONFIGURATION 

The experimental data was obtained using the Air 
Force Research Laboratory low-speed linear cascade 
tunnel (see Fig. 1) at Wright- Patterson Air Force 
Base. The wind tunnel is a modified Aerolab draw- 
down t unnel using a variable frequency motor driving 
an axial flow fan with a diameter of 128.6 cm (50.6 
inches). The tunnel modifications include a new ad- 
justable turbine cascade test section inserted between 
the original tunnel inlet and exit, and a turbulence- 
generating grid. The fan is driven by a 9 3. 2-kilowatt 
(125 HP) motor pulling room air through the 304.8 
cm (120 inch) wide by 266.7 cm (105 inch) high 
bell mouth inlet. The inlet contains honeycomb flow 
straighteners and converges to a test section size of 
121.9 cm (48 inches) wide by 87.6 cm. (34.5 inches) 
high. The tunnel is capable of test section speeds 
from 3 to 20 m/s operating in room air with a turbu- 
lence intensity of approximately 1%. 

To simulate the high turning environment of the 
low pressure turbine, the tunnel features an ad- 
justable test section capable of inlet and exit angles 
of 0 to 65 degrees, independent of each other. For 
the present tests the tunnel is set for the on-design 
specifications for the highly loaded Pratt & Whitney 
PAK-B low-pressure turbine blade with an inlet angle 
of 35 degrees and an exit angle of 60 degrees. As seen 
in Fig. 1, this gives the tunnel a total turning of 95 de- 
grees. The eight scaled-up blades in the test section 
have an axial chord of 17.78 cm (7 inches) with an 
axial chord-to-spacing ratio (solidity) of 1.129. The 
blade spacing is 15.75 cm (6.2 inches) from blade-to- 
blade. Each test blade is molded from solid Ultralloy 
108 white resin and is seven times engine scale in 
the axial direction. To further enhance the approx- 
imately two-dimensional flowfield of the midspan of 
the actual blades, each test blade is stretched along 
the z-axis to 78.63 cm (30.96 inches) high for an as- 
pect ratio of 4.92 to 1. The inner and outer walls are 
shaped to approximate the suction and pressure sides 
of a blade, respectively. There are 8 blades with nine 
passages numbered 1 through 8 starting from the in- 
side bend. Boundary layer measurements were taken 
over the surface of blade 5. 

Measurements of true inlet angle, exit angle, ve- 
locity, and turbulence intensity have been made to 
verify the tunnel operating conditions. The inlet and 
exit angles were measured with a 30 degrees 3-hole 
wedge probe located 1 chord length upstream of the 
blades with a 5.08 cm (2 inch) inclined manometer. 
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All angle measurements were made in the pitch di- 
rection (x-axis) which is perpendicular to the blade 
inlet plane. The measurements showed a close con- 
formity with on-design specifications. The inlet angle 
dropped slightly with an increase in turbulence inten- 
sity. 

Test conditions used in this investigation were for 
low Reynolds number operation at 43,000, 86,000 
and 172,000 based upon inlet velocity and axial 
chord length. These Reynolds numbers correspond 
to 102,000, 205,000 and 409,000 using exit velocity 
and true chord length. The periodicity and unifor- 
mity of the inlet velocity and turbulence intensity 
were measured with a single TSI T1.5 hot wire probe 
suspended from a 3-axis traverse powered by a Dantec 
IMC-2 motor controller approximately 1 axial chord 
length upstream of the blades. Exit velocity and tur- 
bulence were measured with the hot wire approxi- 
mately 1 axial chord lengt h downstream of the blades. 
The wire was inserted into the test section through 
slots denoted in Fig. 1 by S-l and S-4. The 3-axis 
traverse can travel 60 cm (23.62 inches) in each axis, 
providing measurement of the flow properties within 
40.64 cm (16 inches) of the inner wall and 22.86 cm 
(9 inches) of the outer wall in slot 1. Pitch traverse 
movements are measured with a Sony SR50-075A lin- 
ear scale and a Sony LY-51 digital display. Traverse 
movements are accurate to within 0.05 cm (0.0197 
inches) . 

The baseline Pratt and Whitney Pak-B blade was 
tested for freestream turbulence conditions typical 
of high altitude, low speed flight of reconnaissance 
UAV’s. Turbulence intensities between wakes for en- 
gines operating at, these conditions are approximately 
2-4% [4, 5]. Turbulence intensity can be increased 
from 1% to 4% through the use of a square lattice 
grid immersed into the flow 13.17 chord lengths (234 
cm) upstream of test blade 5. The turbulence grid 
consists of 16 vertical and 10 horizontal stainless steel 
tubes, each 2.54 cm (1 inch) in diameter. The tubes 
are spaced 7.6 cm (3 inches) center to center in ac- 
cordance with guidelines suggested by Roach [20] for 
isotropic turbulence. 

Hot wire measurements of velocity are accurate to 
within 0.07 m/s for the high speed case to ±0.24 m/s 
for the lowest speed case. Table 1 displays all of the 
calculated error bands using the method described by 
Kline and McClintok [21]. The error bands are omit- 
ted from the figures for clarity. Turbulence inten- 
sity is determined from the standard deviation from 
the mean of the fluctuating velocity divided by the 
mean velocity. All hot wire data was recorded using 
40,000 samples collected at a sample rate of 10,000 
hz. Temperature measurements are made with J-type 


thermocouples placed at the tunnel inlet, between the 
turbulence grid and the blade cassette, one on blade 
1, and the test section exit. The thermocouples use a 
Kaye Instruments Ice Point Thermocouple Reference 
with the resulting temperatures read into the com- 
puter via the HP 3852A Data Acquisition/Control 
Unit. Thermocouples are labeled T-l through T-4 in 
Fig. 1. 


Re 

U(m/s) 

c P 

0.43 x 10 5 

±0.24 

±0.59 

0.86 x 10 5 

±0.15 

±0.18 

1.72 x 10 5 

±0.70 

±0.04 


Table 1: Error bands for hot wire and pressure mea- 
surements. 

Surface pressures were measured at 40 locations (13 
pressure surface, 27 suction surface) around blades 
4 and 6. The surface pressure tap holes are 0.76 
mm (0.03 inches) in diameter and are located in the 
center 22.86 cm (9 inches) of the blade span. The 
measurements were made using two 48-channel Scani- 
valve J-9 scanners each with a 0.5 psid Druck pres- 
sure transducer. The scanners are powered by two 
CTLR 10P/S2-S6 controllers operated via an HP 
3852 A Data Acquisition /Control Unit and in-house 
software. The transducers are powered by an HP 
6 102 A DC Power Supply and measurements are read 
into the operating computer via the HP 3852A Data 
Acquisition/Control Unit and in-house software. The 
errors in the Cp curves was calculated using the 
method described by Kline and McClintok [21] and 
are listed in Table 1. 

NUMERICAL ANALYSIS 

The flow field is divided into two types of zones 
in the numerical procedure. O-type grids are used 
to resolve the flowfield near the airfoils. The O-grids 
are overlaid on H-grids which are used to Tesolve the 
remainder of the flow field. The thin-layer or full 
Navier- Stokes equations are solved on both the O- 
and H-grids. The governing equations are cast in 
the strong conservation form. A fully implicit , finite- 
difference method is used to advance the solution of 
the governing equations in time. A Newton-Raphson 
subiteration scheme is used to reduce the lineariza- 
tion and factorization errors at each time step. The 
convective terms are evaluated using a third-order- 
accurate upwind-biased Roe scheme. The viscous 
terms are evaluated using second-order accurate cen- 
tral differences and the scheme is second-order accu- 
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RESULTS 


rate in time. Details of the solution procedure and 
boundary conditions are discussed in Ref. [22] 

Two models were used to simulate the effects of 
turbulence. The first model is a two-layer algebraic 
turbulence model based on the work of Baldwin and 
Lomax (BL) [23]. Several modifications were made 
to the original BL model based on previous experi- 
ences with compressor and turbine geometries [19]. 
The second model is a two-equation k — e turbulence 
model based upon the work of Towne el ai [24]. In 
the current, implementation, the equations for the tur- 
bulence kinetic energy and dissipation rate are de- 
coupled from the flow equations and solved using 
an alternating-direction implicit integration scheme. 
Since Newton- Raphson subiterations are used at each 
global time step of the flow solver, decoupling the 
k — e equations from the flow solver should not affect 
the time accuracy of the analysis. The convective 
fluxes in the turbulence equations were discretized 
using first-order accurate upwind differences, while 
the dissipation terms were discretized using second- 
order accurate central differences. The k — e subrou- 
tines were constructed in a modular manner to allow 
the use of different low 7 Reynolds number approxi- 
mations; the Chien low Reynolds number model has 
been used in the current investigation [25]. Previous 
studies have shown that the results predicted using 
the k — e model are sensitive to the specified value 
of the free stream dissipation length scale [19]. The 
dissipation length scale was not determined in the 
experimental portion of this study, so a value w 7 as 
chosen based on experience. 

The low 7 Reynolds number environment in low- 
pressure turbines suggests that the flow 7 may be tran- 
sitional. In the current investigation, natural transi- 
tion is modeled using the Abu-Ghannam and Shaw 7 
model [26]. In the region between the start and 
end of transition the intermittency function is deter- 
mined using the model developed by Dhawan and 
Narasimha [27]. For cases involving separation bub- 
bles the model developed by Roberts [28], and mod- 
ified by Davis et al. [29], is used. 

In the numerical simulations, the O-grid contained 
351 x41 (streamwisex tangential) grid points and the 
H-grid contained 150 x 51 grid points, for a total of 
22,041 points (see Fig. 2). The average value of 
the non-dimensional distance of the first grid point 
above the surface, w r as approximately 0.1. The di- 
mensions of the O-grid were arrived at based on the 
value of and the number of points within the 
boundary layer (approximately 20-25), while the di- 
mensions of the H-grid was determined by performing 
wake convection simulations in the absence of airfoils. 


Three different Reynolds w*ere investigated in this 
study. Re — 43,000, 86, 000 and 172,000. The tur- 
bulence intensity was set at Tu = 1% for all three 
Reynolds numbers. The flow 7 conditions associated 
with the three Reynolds numbers are showrn in Ta- 
ble 2. The exit Mach numbers, the exit flow T an- 
gles and the pressure ratios shown in Table 2 were 
obtained from simulations using the BL and tran- 
sition models. The inlet Mach numbers, inlet flow 7 
angles and Reynolds numbers were similar in the ex- 
periments and simulations. Note, because of the low 7 


Re 

0.43 x 10 s 

0.86 x 10 5 

1.72 x 10 5 

Mi 

0.011 

0.023 

0.046 

Mn 

0.016 

0.035 

0.065 

Qi (deg) 

35.00 

35.00 

35.00 

« 2 (deg) 

-57.92 

-59.03 

-59.26 

Pi/Pn 

0.9998 

0.9991 

0.9970 


Table 2: Flow 7 conditions. 


Mach numbers all the simulations were performed us- 
ing double-precision arithmetic. 

Figures 3 to 5 illustrate time- averaged Mach con- 
tours at the three Reynolds numbers. Figure 3 dis- 
plays a large, unsteady separated flow 7 region on the 
suction surface of the airfoil. The amount of flow sep- 
aration, and the associated unsteadiness, decreases as 
the Reynolds number is increased. There is a long, 
shallow 7 separation bubble in the cove region of the 
pressure surface at all three Reynolds numbers. At 
the low 7 er two Reynolds numbers the form of tran- 
sition was bubble transition, while at the highest 
Reynolds number the transition process varied be- 
tween natural and bubble transition. 

Figures 6 to 8 show time- averaged non-dimensional 
entropy contours for the three Reynolds numbers. 
Similar to the Mach contours, the entropy contours 
show 7 decreasing boundary layer and wake thick- 
ness as the Reynolds number is increased. As the 
amount of separation decreases, the exit flow 7 angle 
approaches the design value of a 2 = 60 degrees (see 
Table 2). The w 7 ake thickness is especially impor- 
tant because it will have a strong influence on the 
performance of downstream blade row 7 s in an engine 
environment. 

Airfoil Loadings 

Figure 9 illust rates a comparison of the experimen- 
tal and predicted blade loadings at Re = 43,000. 
Included in Fig. 9 are results from fully turbulent 
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and transitional simulations with the BL turbulence 
model, a simulation using the BL model with the 
boundary layer tripped at 35% of the axial chord, 
and a transitional simulation using the k — c model. 
The predicted results from the transitional simulation 
with the BL turbulence model exhibit fair agreement 
with the experimental data, although the pressure co- 
efficient is underestimated in the separated flow re- 
gion on the suction surface of the airfoil (implying 
lower velocities and more pressure recovery in the 
simulation). In the experiment the flow separated 
at 60% of the chord and did not reattach, while in 
the BL simulation the flow separated at 63% of the 
chord and did not reattach (see Table 3). The loading 
from the transitional simulation with the k — e model 
shows the same character as the BL solution, but con- 
tains a larger undershoot in the separated flow region. 
Note, however, no attempt was made to tune the dis- 
sipation length scale in the k — e simulations. The 
fully turbulent simulation and the simulation with 
the tripped boundary layer display similar loadings. 
The increased mixing from the turbulent flow sup- 
presses flow T separation and aids in the pressure re- 
covery downstream of the throat. 

A comparison of the experimental and predicted 
blade loadings at Re = 86,000 is shown in Fig. 10. 
Included in Fig. 10 are results from fully turbulent 
and transitional simulations with the BL turbulence 
model, a simulation using the BL model with the 
boundary layer tripped at 35%; of the axial chord, 
and a transitional simulation using the k — e model. 
The results of the transitional simulation with the 
BL model show good agreement with the experimen- 
tal data, exhibiting both the separated flow region 
(beginning at approximately 72% axial chord in the 
experiment and 70%; axial chord in the simulation) 
and the pressure recovery from approximately 87% 
of the chord to the trailing edge. In the experiment 
the flow reattached at 87% of the axial chord, while 
in the simulation the flow did not fully reattach until 
95%) of the chord (see Table 3). The loading from the 
simulation with the k — e model is nearly identical 
to that in the BL simulation up to 75% of the axial 
chord. Whereas the experimental data and BL re- 
sults indicate a pressure plateau from approximately 
75% chord to 87.5% chord, the k — e results show 
pressure recovery from 75% of the chord to the trail- 
ing edge. The fully turbulent and tripped boundary 
layer simulations again produced similar loadings. 

Figure 11 illustrates a comparison of the ex- 
perimental and predicted blade loadings at Re — 
172,000. Included in Fig. 11 are results from fully 
turbulent and transitional simulations with the BL 
model. The results of the transitional simulation 


agree closely with the experimental data. In the ex- 
periment the flow separated at, 75% of the axial chord 
and reattached at 83% of the axial chord. In the sim- 
ulation, the flow separated at 72% of the axial chord 
and reattached at 83% of the axial chord (see Ta- 
ble 3). 


Re 

0.43 x 10 s 

0.86 x lO 5 * 

1.72 x 10 5 

Separation 





0.60 

0.72 

0.75 


0.63 

0.70 

0.72 

Reattach 




(A/fc x ) erpf 

- 

0.87 

0.83 

(A/6 x ) pre<f 

- 

0.95 

0.83 


Table 3: Separation and reattachment data. 


Boundary Layer Profiles 

Velocity profiles were obtained normal to the sur- 
face at 67.2%), 73.0%, 79.3%, 84.8%; and 89.8%; of the 
axial chord. While the predicted results show both 
positive and negative velocities in separated flow re- 
gions, the experimental profiles do not show r negative 
velocities because of the operation of the hot wire. 

The experimental and predicted velocity profiles at 
Re = 43, 000 are shown in Fig. 12. The BL and k — e 
results shown in Fig. 12 include the effects of transi- 
tion. The BL results show good agreement with the 
experimental data at the 67.2%. and 73.0%; locations, 
but the experimental data indicate more rapid growth 
of the separated layer downstream of 73.0%> axial 
chord. The k — e results initially show more rapid 
growth of the separated flow region, but show good 
agreement with the experimental data from 79.3% of 
the chord to the trailing edge. 

Figure 13 illustrates the experimental and pre- 
dicted velocity profiles at Re = 86,000. The BL and 
k — e results show fair agreement with the experimen- 
tal data at 67.2% of the chord, display more rapid 
growth of the boundary layer at 73.0%. and 79.3% of 
the chord, then again show fair agreement with the 
experimental data at 84.8%) and 89.8% of the chord. 

Velocity profiles were not obtained in the experi- 
ments at Re = 172, 000 because of limitations of the 
probe. 

Integral Boundary Layer Quantities 

Figure 14 contains the predicted time-averaged 
skin friction distributions. As the Reynolds num- 
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ber was increased from 43,000 to 172,000 the suction- 
surface separation point (denoted by Cj < 0) moved 
from approximately 63% of the axial chord to 72% 
of the axial chord. The flow did not reattach at 
Re — 43,000, reattaches near the 95% axial chord at 
Re — 86, 000 and reattaches at approximately 83% 
axial chord at Re = 172,000. There is a shallow 
separation region in the cove portion of the pressure 
surface surface at all three Reynolds numbers. 

Figure 15 contains the predicted time- aver aged dis- 
placement thickness distributions. The displacement 
thicknesses on the pressure surface of the airfoil are 
similar for all three Reynolds numbers. The growth 
of the displacement thickness aft of peak suction on 
the suction surface increases rapidly as the Reynolds 
number is decreased (and is caused by increased flow 
separation). 

Cascade Losses 

The experimental and predicted (BL turbulence 
model with transition) cascade losses are shown in 
Fig. 16. Included in Fig. 16 are the baseline losses, 
and the losses with a tripped boundary layer. In the 
experiments, the boundary layer was tripped at 35% 
of the axial chord using a trip wire attached to the 
suction surface of the blade. The location of the trip 
was determined based on requirement that the flow 
remain attached in the region where it was separated 
on the baseline blade. In the simulations the bound- 
ary layer was forced to start the transition porcess at 
35%; of the chord, and the length of transition was 
set equal to 10% of the axial chord. In the baseline 
experiments the losses initially drop off rapidly as the 
Reynolds number is increased aboved Re = 43,000. 
The experimental losses remain nearly constant be- 
tween Re = 86, 000 and Re — 172, 000. The predicted 
losses for the baseline airfoil were greater than the 
experimental values at all three Reynolds numbers. 
The predicted results also indicate a more gradual 
decrease in the loss as the Reynolds number is in- 
creased, which is consistent with the longer separa- 
tion bubbles in the simulations (see Table 3). 

The losses in the experiments with the tripped 
boundary layer are nearly constant across the range 
of Reynolds numbers. The losses are significantly 
lower than the corresponding baseline losses at Re = 
43, 000, but greater then the baseline losses at the 
two higher Reynolds numbers. The predicted results 
for the tripped boundary layer show lower losses than 
the baseline simulations for the two lower Reynolds 
number, and similar losses at the higher Reynolds 
number. These results are consistent with the fact 
that as the Reynolds number is increased the flow ap- 
proaches fully turbulent conditions, and tripping the 


boundary layer should have little effect on the losses. 
The increased losses at the higher Reynolds numbers 
in the experiments may be due to the presence of the 
trip wire. 

CONCLUSIONS 

A series of experiments and numerical simulations 
have been performed for the PAK-B low-pressure tur- 
bine geometry. The trends of the predicted results 
are consistent with the experimental data, and quan- 
titatively the predicted results exhibit fair agreement 
with the experimental data. 

The predicted results and experimental data indi- 
cate: 

• increased losses and flow separation as the 
Reynolds number is decreased 

• bubble transition at lower Reynolds numbers, 
and a mixture of bubble and natural transi- 
tion as the Reynolds number is increased to 
Re = 172,000 

• boundary layer trips are an effective mechanism 
for reducing losses at low Reynolds numbers, but 
not at higher Reynolds numbers 

• in the absence of experimentally determined val- 
ues of the dissipation length scale, the Baldwin- 
Lomax turbulence model (in conjunction with a 
transition model) produces satisfactory results 
for low Reynolds number flows 
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